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In this talk we discuss the calculation of the inclusive semileptonic 
decay lepton spectrum from B -^ XeV. We first discuss the moti- 
vation for renewed interest in this calculation, and discuss the use of 
the OPE and HQET to determine this spectrum. We then compare 
this prediction to that of heavy quark models, in particular that of 
Altarelli, Cabbibo, Corbo, Maiani, and Martinelli IQ (hereafter re- 
ferred to as ACCMM). I will concentrate on i? — > XcSU and not the 
detailed questions of the endpoint of the spectrum which are impor- 
tant for B -^ XuCV. 

There are several reasons for renewed interest in this prediction. 
First of all is the copious production of b quarks at LEP and CLEO. 
A good understanding of the decay spectrum is necessary for under- 
standing fundamental b quark parameters, which are being precisely 
studied, such as the b width and the forward-backward asymmetry 
at LEP and those couplings relevant to the rare decays studied by 
CLEO. Furthermore, with a good understanding of the spectrum one 
can do detailed measurements of heavy quark fragmentation, which 
could give interesting tests of heavy quark theory and QCD predic- 
tions y, ^. Furthermore, in order to do high statistics measurements, 
one needs to do inclusive measurements. On the theoretical side, it 
is of interest to study the inclusive decays because they are under 
much better control theoretically than predictions for exclusive modes, 
where one needs to know hadronic matrix elements. An understanding 
of the lepton spectrum is therefore crucial to extracting fundamental 
parameters. 

In the seminal paper of Chay, Georgi, and Grinstein Q it was 
shown that one can treat the decay spectrum with the operator prod- 
uct expansion and heavy quark methods. They showed the leading 
A/m corrections to the free quark result vanish in the matrix ele- 
ments. Subsequently, Bigi, Blok, Koryakh, Shifman, Uraltsev, and 
Vainshtein P] studied (A/m)^ corrections. They compared the results 
to a free quark, ACCMM model where they did not see the vanish- 



ing of the A/m corrections as manifest. Following, there were further 
calculations of (A/m)^ suppressed effects by Falk, Luke, Savage [^ in 
the context of 6 — > 57 and by Manohar and Wise |p in the context 
of semileptonic decays. We will borrow heavily from these latter two 
papers in our review of the QCD calculation. 

Defining x = -^, e = — j and Xm = 1 — e, the free quark decay 

b 

spectrum is given by 



-^^- = 96vr3 (1 - x)3 [(1 - ^)(^ - 2^) + (^ - "")(^ - ^)] 

(1) 
The idea is then to justify and improve on this free quark result. For 

the study of the inclusive decay, we sum over all final states with given 
quantum numbers. The spectrum and rate of the inclusive decay is 
then governed by short distance physics. Intuitively,the justification 
is that in the ?tt, ^ 00 limit, the decay time is much shorter than 
hadronization time scale. With the OPE and HQET, this can be made 
precise. The crucial observation is that there is a large range of mass 
scales in the final state m'jj < P^ < m^ so that the energy flowing 
through the hadron system scales with the mass of the decaying heavy 
quark. In the TTi ^ 00 limit, this is much larger than the QCD scale. 
Away from P^ ~ ^-i)) the internal quark is far from mass shell. In 
this case, we expect the OPE to be useful; it will prove valid except 
near the endpoint. 

The general idea is to relate the square of the matrix element of 
interest to the imaginary part of the time ordered product of currents. 
One can then perturbatively compute the time ordered product with 
the operator product expansion (in the region of phase space where 
it is valid). However unlike a standard OPE it is useful to expand 
in the heavy quark mass rather than q^ in order to apply HQET. 
The coefficients of the heavy quark operators are then obtained by 
evaluating the time ordered product between quark and gluon states. 
The time ordered product is evaluated in the end by taking matrix 
elements of the operators appearing in the OPE between meson states. 



One can use the HQET to learn about these matrix elements. 

Semileptonic decay is determined from the weak hamiltonian den- 
sity: 

Hw = -Vjb ^ qjj^'PLb ej^Pii^e = -Vjb ^ Jpl^., (2) 

The inclusive differential decay rate for H^y -^ Xu^c^T^ (here we will 
be restricting attention to meson decay) is governed by 

I^f = {27rfY.^' (m -Q-Px) {H,{v,s)\jf \X) {X\ JJ \H,{v,s)) 

X 

(3) 
where 

PH, = Muy = mbv^" + k^ (4) 

The W'^'^ can be expanded in terms of five form factors 

W^"" = -g'"'Wi+v^'v''W2-ie^"''"^Vaq(}W3+q''q''W4+{q^'v'' + q^v^) W5, 

(5) 
One then relates this to a time ordered product via Im T^'^ = 
—ttW^'^ where 

•' s 

= -g'"'Ti + v^v''T2 - ie'"''''^Vaqi3T3 + q^q^T^ + {q^v" + q^v^) Tg 

Let us consider the analytic structure. There are cuts correspond- 
ing to the decay of interest, to the process eVH -^ X, and to e^i^eH — > 
X. The idea is then to do a perturbative calculation away from the 
physical cut, which is always possible away from the endpoint. 

As an example, we find the leading order result. This is done ex- 
plicitly in [y|. The matrix element of the time Fourier transformed, 
time ordered product of currents 



IS 



{mfjV — q + kY — rrij + ie 



-i /"d^xe-*«-^r(J^tj^) (6) 



n7^PL {mbii-i + ^ + mj) j" Pl u, (7) 



Using 

^t^^^Y = g'^'^j" + g'^'^-f'' - g^^'i'^ + ie^"''"^jp-/5- (8) 

one finds the order fc" term is 

-— u I {rubV - qY Y + {iTibV - qY 7^ - {mi,i) - f) g^" 

-ie^^'^-P{m,v-q)^^p]PLU, (9) 

where 

Ao = {rrihV — qY - rrij + ie (10) 

By replacing the matrix element with the b quark operator which 
yields this amplitude when evaluated between b quark states, one ob- 
tains 

-^ {{m,v - qYg"^ + {m,v - qY g^"" - {m^v - qf g^' 
Ao >- 

-ie^'-^^{m,v-q)^}b^xPLh. 
Because it is a current, one has 

{Hb{v,s)\b^^h\Hi,{v,s)) = v^, (11) 

since b quark number is an exact symmetry. 

So we can evaluate the leading order contribution exactly. By 
taking the matrix element of the OPE between H states, we get 

TO = i-m, 

So to get the amplitude W^^^ which determines the spectrum, we need 
the imaginary part of T^^^ . This is readily obtained from 

6{imbV - qY - m'^) (12) 

-6'{{mi,v-qY -m^) (13) 



Ao 
1 

^0 



i_ = ^s"{{m,v-qY-m') (14) 



^0 



(Here we only use the first one for the leading order term) 

The mass shell condition above when integrated over phase space 
leads to just the free quark decay distribution! 

However, with these methods, one can also consider higher order 
corrections. The important result will be that Aqcd/"^ corrections 
vanish and only two matrix elements, one of which is known, are 
required to obtain (Aqcd/'"T')^ corrections 

We first introduce heavy quark fields and the heavy quark effective 
theory. Recall that in the heavy quark limit, the heavy quark mass 
and spin decouple from the soft degrees of freedom. The heavy quark 
state looks like a static point source of charge. The b quark is almost 
on mass shell p = nibV + k. 

Define the heavy quark field 



[1 + ^) 



^inibV-x^ _|_ 



(15) 



with leading order lagrangian 



C = b^ {iv • D)by + 



(16) 



Now expand in terms of HQET operators. We follow [0] who give 
the leading term in terms of the standard quark field, but the mass 
suppressed operators in terms of heavy quark fields. 



rub 



Zmh 



4ml 



O2 + 



(17) 



O2 



bTb, 

by T iD^6„ , 

by T iD^iD,^6„ 



As before, at leading order, we have 



(18) 
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Now consider the first order mass suppressed terms. 

{M\Oi\M) = {M\hTSD^h\M) = {M\hTv^v ■ 'iBh\M) (19) 

This vanishes by the equation of motion! Notice this critically depends 
on our phase choice in defining heavy quark field. 
An arbitrary counterterm 

6mKK (20) 

was taken to vanish. With this choice of quark mass, all linear correc- 
tions to the free quark result vanish! 

At higher order in l/m, one needs the two matrix elements 

Kb = -{Hb{v,s)\by ^ by\Hb{v,s)) 

Gb ^ Zb{Hb{v,s)\b J^f'f K\Hb{v,s)), (21) 

The second matrix element is known because it is related to the 
known spin dependent mass splittings. The first is not known, and is 
a parameter to be determined. 

So, to summarize, the leading order result for the rate and spec- 
trum is the free quark result. The first order corrections vanish. At 
second order, there is one unknown coefficient and one which has 
been determined. This analysis has assumed we are far from end- 
point, where the OPE is valid. Had we applied this near the endpoint, 
higher dimension operators would not be suppressed and it would be 
necessary to average the spectrum over a range of energies. 

One can compare this result to the predictions of models. It is 
interesting to see how models reflect and differ from these general 
QCD predictions. Consider for example the ACCMM model. Here 
I review work of Csaba Csaki and myself [P], but other interesting 



references are that by Grant Baillie [^ and Ref. [10|. One finds that 
one can always define a quark mass so that l/m corrections vanish. 
However, at 1/m^, the model would differ from the QCD result. 
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In the ACCMM model, one models a B meson decay as disinte- 
gration. There is a spectator quark with mass rrigp and momentum 
distribution 0(|p|). 

The b quark momentum is determined by kinematic constraints. 
The lepton spectrum is determined from the decaying off shell b quark. 
Define 



Ew = Mb- P^ + ra% (23) 

The invariant mass of the b quark will then be 



W^^ = m| + vr? - 2mB \ P^ + m'j (24) 



Define 



2Ee 

X = 

m 

2 

mi 
mf 

Recall the formula for free quark decay. 

■v2 ^5 ^2f^ ^\2 



d. - 96.3 n_,.3 [(l-x)(3-2.) + (l-x„)(3-.)] 



dT{nn„E) _ G^pw? x\xm-xY 

(25) 
Then the lepton spectrum in the ACCMM model is 

dcosO 



dE Jo -^1-1/ _/ ^ dE'dcosO 

f d cos 9„ ^, ^ ^. ^. 

/ — :^5{E - lE' - j(3E' cos / 



dTB _ [ , 2^n u 1 /■ ^r , AE' 



dE -J ^^^^^(H)wi ^iE^d^e'^'^'^ ^''^ 



Using j3^^ = pEw /W{py yields 

-^ = J dpp ^i\p\)7;z^r: L -^ — TE^, — (27) 

where 



dE J "^"^ "^'"^"IpEw Je^ E' dE 



EW 



E^ = 


Ew^P 




(28) 


-t^max - 


W 
= min{E+,—{l- 


- Xm)} 


(29) 


Pmax - 


TUb ITT-} 

2 ml 




(30) 



Now do a heavy meson mass expansion. 

-dE=Jo ^^^'^^' ^i^E^^S^ Je -E^[w) ['-w) 

(31) 

We can exphcitly evaluate the E' integral, to get 

(32) 
Take pf small. 

dTs _dr% G'j.E'^Ml ( E \p/ j:>(PF.2d^q 

1[e~1[e^^[2^\/^V^~)^^^^^ He ^^^^ 

It looks like there are nonvanishing linear corrections, but that is be- 
cause we expanded in the meson mass tub , rather than the quark mass 
rrH). To do a heavy quark expansion, one needs to define a quark mass. 
This can be done so that the linear correction vanishes. 

dVB _GlE^ fP™^ dp ( 2 opp2 8 2pA 2,., n 

(34) 
If we define 

rrib = {Wip)) (35) 

Up to quadratic terms in p,we have 

(W) = {Ew) = nib = TUB- {p) (36) 
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{f{W,Ew)) = f{mB,mB)-{p)f'{mB,mB)+0{p^) = f{mh,mb)+0{p^) 

(37) 

We see that it was possible to define a b quark mass so that Hnear 

corrections vanished because the p dependence was all through Ew. 

So we can eliminate linear corrections through a proper choice of mb- 

It is easy to extend this to nonzero spectator mass. 

(Ew) =mB- {^m% + p^) (38) 

Now consider the total rate and the full spectrum. 

(39) 
However, the spectra themselves are very different near the endpoint. 



dE Jo '"^ ^^"-"2p Je_ E' 

r"""^^ 2^/1 ,.W /-^i dT{W,E') 

2 

where Ei = W{1 - e)/2 , pi = mij/2(l - ^ - £; for m^p = 0. This 

rag f 

spectrum extends up to E^ax = "'""a"'" [} ~ (mB-4,p)^ 

We see the ACCMM spectrum and the free quark spectrum deviate 
due to two effects near the endpoint. The first term is not integrated 
up to p ~ oo and the second term integrates over a rapidly falling 
spectrum so it is no longer a good approximation to replace {fiW)) 
by /((^))- The spectrum extends beyond naive quark mass endpoint. 
Since we have a specific model, we can also investigate the question 
of how large an energy interval must be averaged over to get good 
agreement between the model and the free quark prediction. In ref. 
[0] we used the averaging function 

^r, , f 1 -{^-£^^)dT{E), , , , 
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We expect to require an average of approximately 2pf to get (pf/m)^ 
agreement. This was appproximately correct. We also checked that 
the best fit quark mass is generally very close to the exact quark 
mass, so the deviation of the model from the free quark prediction is 
probably not very important for practical purposes. 

At higher order, the predictions of this model will not agree with 
general QCD predictions. This is because gauge invariance is not 
incorporated. There are corrections to the free quark result due to 
the fact that terms proportional to (pq) and {ff) both yield l/m^ 
corrections. Recall that in the heavy quark theory the operator (Dq) 
contributions only at higher order in the heavy quark mass expansion 
by the equation of motion. In the ACCMM model, both contribute 
at the same order. Obviously, the ACCMM model also doesn't give 
you spin dependent gluon operator. We conclude that although the 
model is probably adequate in practice (as would be free quark decay) 
it does not properly incorporate QCD dynamics. 

We conclude that there has been much advancement in our under- 
standing of the lepton decay spectrum. It is unclear at what level this 
will be tested. Most of the large deviations from the free quark predic- 
tion are in the endpoint region, where one must smear to get reliable 
predictions. Furthermore, the higher order correction involves an un- 
known parameter and requires very accurate measurements. However, 
from a practical point of view, when using the spectrum to extract b 
quark couplings, we see that the free quark decay models b quark de- 
cay very well. It might be better to fit the spectrum to a free quark 
spectrum that that of a model, such as the ACCMM model. The 
vanishing of Aqcd/'"T' corrections is very important to this conclusion. 

Acknowledgements: I thank my collaborator C. Csaki. I also 
thank the Institute for Theoretical Physics in Santa Barbara where 
this writeup was completed. 
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